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We investigate the twist two gluon distribution functions for spin 1/2 hadrons, emphasizing 

^»0 , intrinsic transverse momentum of the gluons. These functions are relevant in leading order in 

the inverse hard scale in scattering processes such as inclusive leptoproduction or Drell-Yan 

scattering, or more general in hard processes in which at least two hadrons are involved. They 

show up in azimuthal asymmetries. For future estimates of such observables, we discuss specific 

bounds on these functions. 
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' ^ ■ Gluon distribution and fragmentation functions are fundamental quantities in the study of deep inelastic 

I ' scattering processes. In fact, together with their quark and antiquark analogues, these process-independent 

M-i, quantities describe the soft parts of the scattering or, in other words, the deep structure of the hadrons. The 

pi ' partonic (distribution and fragmentation) functions cannot yet be calculated from first principles because we 

lack the non-perturbative treatment of the strong interactions. However, valuable information on these functions 

can be obtained via lattice calculations or theoretical models. 

As soon as more than one hadron is involved in a hard scattering process, it is essential to take into account 
the transverse momentum of the partons. For instance, transverse momentum dependent quark distribution and 
C^ ' fragmentation functions show up explicitly in several semi-inclusive cross sections, in particular in azimuthal 
asymmetries. In the calculation of QCD corrections for these cross sections, the inclusion of transverse mo- 
mentum dependent gluon distributions and fragmentation functions will be necessary. This is our motivation 
to study in this paper the transverse momentum dependent gluon functions. We will follow the correspond- 
ing treatment for quarks developed by Mulders and Tangerman |l|,^, following earlier work by Ralston and 
Soper@. 

In a diagrammatic expansion of the hadronic tensor in powers of the strong coupling, one finds an infinite 
number of gluon correlators, which are essentially matrix elements of non-local products of gluons fields (and 
sometimes quark fields) between hadronic states. The simplest of these matrix elements are the ones that 
contain only two gluon fields. In the ^4+ = gauge, they completely define the twist two functions through the 
appropriate choice of Lorentz indices and projections. One has to make sure that the starting point is a gauge 
invariant object, which turns out to be a non-trivial matter. 

More complicated correlators, namely with three gluon fields, must also be studied. Some of these correlators 
will precisely provide the link operator needed to define gauge invariant functions, and others will reduce to 
gluon-gluon and gluon-quark correlators using the QCD equations of motion. We will discard all correlators 
that contribute at order 1/Q^ or higher in the cross sections, Q being the hard scale. 

The paper is organized as follows. In the next section we start from the gauge invariant gluon-gluon correlation 
function and derive all twist two and three gluon distributions and we make an analysis of the link operator. 
We show that the link appears naturally when a certain class of diagrams is summed leading to a specific path. 



In section III we introduce a specific helicity basis for nucleons and give the spin representation of the twist 



two part of the gluon-gluon correlator, which can be used to derive the natural interpretation of the lower twist 
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FIG. 1. The giuon correlator 



functions as giuon densities in the framework of the parton model. In section IV the formalism is extended to 
the giuon fragmentation functions, after which we discuss some bounds. We end up with some conclusions and 
suggestions for future investigations that can use the formalism developed in this paper. 

II. GLUON CORRELATION FUNCTIONS 

In order to connect gluons in a hard scattering process to hadrons appearing in the initial or final state, we 
will use (lightcone) correlation functions m-vl ■ Our starting point is the correlation function 

S'^''{k;P,S;n) = / -^ e*«(P, 5|A''(0)A^(e)|P, 5), (1) 

J (Stt) 

diagrammatically represented in Fig. w. The vectors P and S are respectively the momentum and the spin of 
the hadron, while k stands for the momentum of the giuon. Additional dependence can come for instance from 
fixing the gauge using a vector n. A summation over color indices is understood. If, as done in this paper we 
use the notation A^(C) = A"(^)r°, where the T" are the generators of the SU{3) color group, this summation 
is actually an appropriate tracing. 

A corresponding gauge invariant object is the quantity 

r^'^--P-{k;P,S) = J 0^e'^< {P,S\F^''iO)U{0,0 FP^0\P,S), (2) 

where Fp,y(^) = F'^^,{^)T'^ is the field tensor, related to the potential by F^^ = d^A^, ~ d,jA^ — ig [A^^, A^,]. The 
link operator hi will be studied in detail in one of the next sections. 

A. The Lorentz structure of the giuon correlator 

The Lorentz structure of the giuon correlator is limited by constraints following from hcrmiticity and parity 
conservation. These are 

TP''■'^"'*{k■P,S)^T^''''P''{k■,P,S), (3) 

where k = {k'-\ — fc'). A possible parameterization of F^'^-p'^ compatible with these constraints is 
TP'^-^P'^ (k; P, S) ^ Xie^"""^ e"" c.b + ^P^'' g"^^" P''^ + ^ k^^" g"^^" k"^ 

M2 ^ M2 ^ 
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+*?FT e^''P^''PP\k-S) + i^ e^'"'''^'' kP\k ■ S) 
+1^ et^'^P^'^ kP\k ■ S)+i^ eM-fekpp](fc . s) 

M3 ' M3 

X20 + 1X21 ^^^fep p[p_5.o-] I ^20 - ^^21 ^pcrfcP plp^H 

M3 M3 

X22_f_iX23 ^^^fep ^[p^ff] X22 - iX2y. ^p^^P ^[t^ gu] 

M3 Af3 

+ ^''m5^'' ^p.feP fc[pp.](fc . ^) ^ -^24-^»X25 ^p.,P fc[MpH(fc . S). (5) 

The amplitudes X^ as well as the original correlator T'^'^'P"' {k; P, S) have dimensions Af^^. With the chosen 
parametrization, i.e. the appropriate introduction of factors 1 or i for symmetric and antisymmetric tensors, 
hermiticity in Eq. (pf) implies that one finds real amplitudes Xi. The parity constraint in Eq. (m requires 
that even numbers of e-tensors are combined only with vectors k and P, while odd numbers of e-tensors are 
combined with the axial vector S besides vectors fc, P. Since S parametrizes the nucleon density matrix, it can 
only appear linearly. 

Time reversal invariance, when applicable, imposes a third condition, 

rp,,p^(k;P,S)=T^''--''"*{k-P,S), (6) 

which implies X* = —Xi for the amplitudes X5, Xy, Xiq, Xis, X20, X22 and X24. For this reason, such 
amplitudes are called T-odd. They thus vanish when time-reversal can be used as a constraint. 

B. Twist expansion 

Leading and non-leading contributions to the hadronic tensor are easier to identify if one uses a suitable 
parameterization of the hadron momentum and spin vectors in terms of two light-like directions, n+ and n_ 
(such that n^-n_ = 1), and two transverse directions. They are chosen such that the hadrons have no transverse 
momentum, what means that P can be written in terms of the light-like vectors. The momentum of the gluon 
and the spin vector of the hadron must include a transverse component: 

A'P 
P = P+n^ + ^j^n., (7) 

p + fe^ 
k = xP+n+ + - „ 7 n_ -I- fc^, (8) 

P+ M 

S ^Sr^— n+-Sr^^^ n_+ST. (9) 

The quantity x represents the fraction of the light-cone momentum in the -I- direction carried by the parton. 
The parameter S^ and the two component vector St are such that 5*^ -I- S^ = 1. The quantity S^ is referred 
to as the helicity. Having defined the n± vectors one has transverse tensors g!^'^ and e^'^ defined as 

5^'' = g'--<<-<n^, (10) 

et^" = £"+"-'"' = e'+f"". (11) 

The expansion in lightlike vectors shows its usefulness only when the correlation functions are used in a calcu- 
lation of a hard scattering process in which a hard timelike or spacelike vector appears setting the hard scale. 
An example is the momentum transfer squared, q^ = —Q^ in deep inelastic leptoproduction. For the soft part 
in the process, described with the correlation functions it implies that after all calculations are finished, the 
P+ is of the same order of magnitude as the hard scattering scale Q. Simple power counting tells us that the 
most important contributions from F'^''''"^ are the ones with the largest possible number of + indices. Due to 



the antisymmetric character of the field tensor, these are obviously r+''+-' (referred to as twist two), followed 
by r+*'^ and r'^''+ (referred to as twist three), where i,j,l, ■ ■ ■ indicate transverse indices. 

For deep inelastic scattering processes we always need the soft parts integrated over the momentum component 
k~ . Starting with the twist two part, we define 



Mr^(j;,fcr)= I dk- V+''+^{k;P,S)= f 
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This quantity depends on the momentum fraction x and the transverse momentum k^ besides the (suppressed) 
dependence on the target momentum and spin, and is conveniently expressed in terms of transverse tensors 
and vectors. Concerning the dependence on the hadron spin, we furthermore separate the unpolarized (O), 
longitudinally polarized (L) and transverse polarized (T) situations. This leads to 



i q{x^ krp) 



X P+ 
2 "m 



-g^G{x,kl) 



^9'4 



V 



Krj 



ry(x,k^) = - ^- 



X P+ 

2 "m 



-ie^S^ AGLix,kl) + - 



h.i h.3 

'm^ ' ^^ 2M2 



H^{x,kl) 



2M2 



S,AHt{x,k'^) 



(13) 
(14) 



J- rp \Jj. t\j rj- J 



X P+ 

2 Id 



^UtS^ 



I] ^j_ 



M 



Gxixjk'^) — ie^ 



K'j^ ' J J 



M 



AGrix^kl) 






+ - 



2Af2 M 



AM 



AHrix.kl) 



2M- 



■AH^{x,kl 



(15) 



where the expressions of the functions in terms of the amplitudes Xi can be found in the appendix. The factors 
in this parametrization are chosen in order that G can be interpreted as the giuon momentum density |q] , which 
will become clear when we discuss sum rules at the end of this section and in the next section. Actually also the 
use of the combination AH^ = AHt - AH;^'-^\ where AH^'-''^ = {kl/2]VP)"'AH^ (with similar definitions 
for other functions), is done because it is nicer for interpreting the functions. 

When we use the soft parts in calculations up to 0{1/Q) we need r+*'^ and r*''''+, again integrated over 
dk^ , which we refer to as twist three contributions, 



Mr-(x,fcr)= f dk- r+*'+^(fc;P,S'), 
Mr^-\x,kT)= [ dk- r^''+(fc;F,S'). 
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They are again parametrized in terms of a number of functions. We obtain for the various hadron polarizations. 
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Again the functions expressed in terms of the amplitudes are given in the appendix. While the functions for 
twist two are real functions, those for twist three are arranged in terms of complex functions in such a way that 
the T-even functions correspond to the real parts and the T-odd functions correspond to the imaginary parts. 
For twist two, the functions Gt, ^H^, AH^, and AHt are T-odd. 

Let us make a short remark on the names of the functions, which follow for the indices in part the notations of 
quark distributions as introduced in g (and extended in [Qj^). The distributions are represented by G, H, AG 
or AH. The names G and AG are reserved for functions that do not involve uncontracted momentum indices. 
These do not flip the gluon helicity, and represent unpolarized (G) and polarized (AG) gluons, respectively. 
The functions H and AH flip gluon helicity in unpolarized or polarized targets respectively, as we will discuss 
in detail below. We distinguish between the longitudinally polarized spin 1/2 target being multiplied by S^, 
which acquire a subscript L and those appearing in a transversely polarized spin 1/2 target being multiplied 
by the transverse spin of the hadron, which acquire a subscript T. If there is an uncontracted component of 
the transverse momentum of the gluon multiplying the function or if needed to avoid double names, we add a 
superscript _L. Finally, twist three functions are given an additional subscript '3'. 

In the next step we perform the integration over the transverse momentum of the gluon to arrive at the 
distribution functions which depend only on x and are important in deep inelastic inclusive measurements. We 
define 



r'{x) = j £k^ r'{x,k^) = y ^ e^'=-« {P^S\F+\Q)F+\0\P,S) 



(24) 



and similarly for F* {x) and F'^''(a;). 
We find (combining the polarizations) 



X 



r^^(x) = I ^ [-g^G{x) - S,i^ AG{x)\ , (25) 

r-(x) = I *e^.^ AG3t(x), (26) 

V'^\x) = I i4s'^ AH^t{x). (27) 

where G{x) = J d^k^ G{x,kl) and similarly for AGst and AH3T, while AG{x) = ^ (fk.rAGL{x,kl). The 
functions AG3T and AH^t are in essence the functions Hi and H2 of Ref. ||] . 

C. Sum rules 

Local hadronic matrix elements are obtained from the gluon correlation functions after integration over dk^ , 
e.g. 

^j2py ^ j^j I ^^+ p»j(^^ ^ i^p^ S'|F+^(0)F+'^(0)|P, S). (28) 

The trace of this quantity is precisely the gluon part of the energy momentum tensor. Using the parametriza- 
tion of F'^ (x) in terms of the gluon distribution G{x) one finds 



M^F" = 2{P+f [ dx xGix) = {P,S\e++\P,S). 
Jo 



(29) 



To derive this last relation we used the fact that the integral over dx has a support between —1 and +1 and 
the symmetry relation G(— .x) = —G{x), which follows from the commutation relations for gluonic fields. The 
number 

eg= f dx xG{x) (30) 

Jo 

thus is identified with the fraction of light-cone momentum carried by the gluons, < e^ < 1. 
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FIG. 2. Some contributions to the fictitious highly virtual ferniion-nucleon cross section 

D. Gauge Invar iance 

Of course the object appearing in the diagrammatic expansion, S^^{k; P, S), is not gauge invariant. We solved 
half of the problem by starting with r^"-'"^. In particular, in the A+ — gauge one has F+* — d'^A!^ and thus 
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{k;P,S)^-{k+f S'^k;P,S). 



(31) 



In a general gauge, however, one also needs to consider matrix elements of the form (AM+A-'), (A^A'^A'^A^), 
etc. Two simple leading contributions are shown in a fictitious highly virtual fermion-nucleon scattering process 
in Fig. 0. These will contribute at the same order in an expansion in the inverse hard scale. They will assure 
that in a general gauge one also finds other terms such as [A"*", A'lterms in F+% and more importantly a gauge 
link operator. To be precise in a fictitious calculation as in Fig. eI one finds that the field A^(^) appearing in 
the correlator in Eq. [| is to be replaced by 
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where 77+ = ^+ = 0, and tjt — £,t and 



U{oo,Cj ^V exp -ig / ds^'A^{s) 



(33) 



where the path runs along the minus direction from the point ^~ = cx) to ^~ with ^+ = and (^^ fixed, 
analogous to the path for quark correlators 19] . In the fc^ integrated correlators both links run along ^^^ = and 
a straight link between the lightlike separated points in the (F+'(0) F+-'(^)) correlator in Eq. 24 remains. For 
the non-integrated correlator in Eq. O the links do not close, but with the physical assumption that hadronic 
matrix elements of the type (A^ (0)^^(77" = oo)A^(^~)) vanish this does not pose a problem. Furthermore 
when considering weighted fc^^-integrated cross sections as was for instance done for quark field correlators in 
Ref. pOJ, one anyway reduces the matrix elements to lightlike separations. 

We note that in different processes, e.g. lepton-hadron scattering or Drell-Yan scattering the paths will run 
from different points, that is ^~ = 00 and —00, respectively. In A^ = gauge, the difference, however, is 
precisely A^(^~ = 00) -t-A^(^^ — —00), a quantity remaining to be fixed to fully fix the lightcone gauge [OJ12I. 



III. THE TWIST TWO FUNCTIONS AS DENSITIES 



The fact that we in the gauge A'^ = are left with S"^, a matrix element bilinear in the gluonic fields, suggests 
that it might be possible to find a probabilistic interpretation for some distribution functions. This has been 
discussed in detail in several papers, expanding the transverse gluon fields in modes. We follow here a slightly 
different route that allows us to draw conclusions on the newly introduced leading twist transverse momentum 
dependent functions of the previous section. 

The basic idea is the observation that M*-' = (2M/a;P"'")r*-' is a two by two matrix in the two transverse 
polarizations that for any diagonal element is a (positive-definite) density. Generalizing also to a matrix in the 
hadron spin space, one has 



= 5](P„|F+X0)|P, A)*(P„|F+^(0)|P, A') S (P+ - (1 - x)P+) . (34) 

n 

In principle it does not matter for our considerations if we use the matrix for (real) linear polarizations or 
for the circular polarizations of the gluons. For interpretational purposes, the latter however is more common. 
Using the circular polarizations 

\±)=T^{\x)±i\y)), (35) 



we obtain the matrix elements 



M++ = i(Af" + M22) - ImM^^, 
M+- = --(M" -M22)+i7^eM^^ 
M-+ = -i(M" - Af22) - i7^eMl^ 
M— = -(M"+M22)+ItoM1^ 



Explicitly we find from the parameterization in Eqs |l3| - |l5| the matrix elements (in gluon polarization space), 
M++ = G + Sr, AGl - S^ %^(sin(/. Gt - cos0 AGt) + S^ %^(cos(/) Gt + sin(/) AGt), (36) 
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-«e ' - -r^ AHt{S^ + iS^), (38) 



M— = G-S^ AGl-S^ ^(sin0 Gt + cos AGt) + 5'| ^(cos(/) Gt - sin (/. AGt). (39) 

In order to make the nucleon spin explicit we use the connection 

r^(x, K;P,s)^J2 ^a'a(s) rAA'(^: fe-; ^: ^)' (40) 

A,A' 

where p(S) is the spin density matrix for a spin 1/2 particle characterized by the spin vector S = (5*1,, St), in 
its rest frame given by 

p(S) = i(l + S.<T). (41) 



Using the explicit form of the density matrix we can write (Mfl) as 

r«. = i (r^+ + ri_) + ^ (r^:^+ - r^_) + ^ (r^ + ri+) + ^ [rX^ - r%) . (42) 

We shall now apply this to the matrix AP-', where we include also T-odd functions. We then obtain a 4 x 4 
matrix in the gluon (g) nucleon spin space [\gluon; nucleon) basis | + ; +), | + ; — ), | — ; +) and | — ; — )], 
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FIG. 3. The gluon correlator 
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(43) 



The matrix representation is also convenient to find the physical meaning of the distributions. Well known is G 
which measures the number of gluons with momentum (a;, fey) in a hadron. The functions AGl (AGt) represents 
the difference of the numbers of gluons with opposite circular polarizations in a longitudinally (transversely) 
polarized nucleon. The off-diagonal function H-^ also is a difference of densities, but in this case of linearly 
polarized gluons in an unpolarized hadron. Using the circular polarizations, H^ flips the polarization. 



IV. GLUON FRAGMENTATION FUNCTIONS 



The procedure to analyse the gluon fragmentation functions is very similar to what has been done for the 
distributions in the previous sections. 
Using as parameterization of the vectors 
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for a gluon with momentum k fragmenting into a hadron with momentum P^ and spin S^, we consider the soft 
part (see Fig. S) 



5^'^(fc; P,, S,,) = ^y A^ e^^-« {0\A''iO\Ph, S^-X)(P,,, S,,;X 
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or the appropriate gauge-invariant object 



f^''^^-(fc;P,.,5,.) = ^y'^e*«(0|F''-(O|P.,5„;X)(P;„5„;X|Z^(C,0)F^^(0)|0). (48) 

For the description of fragmentation in leading order in the inverse hard scale, we need this correlation 
function integrated over one lightcone direction, with the above choice for P/i, being the momentum k~^ , 

(49) 



M^r'J(z,fc^) = / dk+ T-r--\k-Ph,Sh). 



which, separating the polarizations, is parameterized as 
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where the argument of the fragmentation functions, k'^ = —zkj. is the transverse momentum of the hadron 
with respect to the gluon. The factors are chosen such that JdztPk'^ zG{z,k'^) = {z)^ is the fraction of 
momentum of the struck gluon taken by the hadron h, for which we have '^i^{z)h = 1- While for distribution 
functions T-odd functions might appear via special mechanisms dealing with initial state interactions or gluonic 
poles, this is not the case for fragmentation functions |13|Jl4[| , where time-reversal symmetry cannot be used 
as a constraint because of the explicit appearance of out states \Ph]X) in the definition. Thus one expects 
nonvanishing fragmentation functions Gxizjk'^), AH^{z,k'^), AH^{z,k'^) and AHT{z,k'^). 



V. BOUNDS ON THE DISTRIBUTION AND FRAGMENTATION FUNCTIONS 



As discussed before, we can organize the distribution functions in a matrix representation in the gluon 
nucleon spin space. For distribution functions one has (omitting the T-odd functions) a matrix 
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Requiring any diagonal element to be positive gives using the diagonal elements the trivial bound 

\AGL(x,kl)\<G{x,kl). 
Using all possible 2x2 submatrices, positivity leads to bounds 
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These bounds can still be sharpened by using the eigenvalues of the full 4x4 matrix in analogy to what was 
done for quark distributions in Ref. |15| . 

For fragmentation functions the matrix {Mh/ P^) r(2:, fc^) contains the various fragmentation functions, now 
including the T-odd functions. It thus is the same as the matrix in Es, but with hat functions depending on z 
and k'. The constraints become 
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FIG. 4. Gluon distribution and fragmentation functions contributing to the single spin asymmetry in p^p -^ ttX 
scattering cross section. 
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VI. CONCLUSIONS 

We have given a full classification of the gluon distributions and fragmentation functions relevant in hard 
scattering processes in leading order in the inverse hard scale including transverse momentum dependence. 
Some results for subleading (twist three) correlation functions have been given also. The inclusion of trans- 
verse momentum dependence is needed in processes involving at least two hadrons. Examples of such processes 
are 1-particle inclusive leptoproduction or Drell-Yan scattering. In these processes one can become sensi- 
tive to transverse momentum dependence, in particular when one considers azimuthal dependence in the final 
state |p],p| jlO|jl^ . We note that in electroweak processes, the gluon correlation functions do not enter at tree 
level but only at higher order in a^ ■ This also implies their relevance in the study of evolution of the transverse 
momentum dependent flavor-singlet quark distribution and fragmentation functions. Also in related processes, 
e.g. X2-production in hadron-hadron scattering the relevance of gluon distribution functions has been empha- 
sized and investigated [|6| . Since we have discussed both distribution and fragmentation functions we have also 
classified the T-odd functions, important in the latter case. In single spin asymmetries at least one of the 
functions describing the soft physics is a T-odd function. 

In order to illustrate the importance of also considering gluons, we consider the contributions to pion pro- 
duction in pp^ scattering in which one of the protons is polarized ||l7|-pi||. In this process a large single spin 
asymmetry is found. We note that gluon correlation functions can play an important role here. This is another 
way of looking to the approach in Rcfs. |lq-[2l|]. Since in this case several mechanisms are considered in these 
papers, we need both T-odd and T-even distribution and fragmentation functions. In Fig. H three contribu- 
tions are shown of gluon correlation functions producing an asymmetry. Actually kinematics requires one of 
the three partons to be off-shell, which means that we need the asymptotic transverse momentum dependence 
of the soft parts, i.e. the evolution of the soft parts. Nevertheless the structure of the gluonic soft parts is 
sufficient to indicate that diagram (4a) will produce a sin((/)^ -I- ips) asymmetry (being equivalent to the Collins 
asymmetry in leptoproduction). This asymmetry is proportional to {hi{xi) x H^{x2)) x h^{z) involving the 
transverse momentum dependent gluon distribution function H^, the transverse spin distribution /i° and the 
Collins function h^ (we use here the hat-notation for the fragmentation function in order to avoid confusion 
between fragmentation functions and gluon distribution functions). Diagram (4b) will lead to a sin(07r — </>s) 
asymmetry proportional to {Gxixi) x /f (X2)) x Df{z) involving a T-odd gluon distribution function, the unpo- 
larized quark distribution function /f and fragmentation function Df. Diagram (4c) gives a similar asymmetry 
as (4b) with the unpolarized gluon distribution function G and gluon fragmentation function G. 
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APPENDIX A: GLUON DISTRIBUTIONS 

Leading and subleading gluon correlations are distinguished via the Lorentz indices. It is therefore convenient 
to rewrite the covariant expression in Eq. p^in terms of {n^, n_,kT., St)- We start with the correlation functions 
with the maximal (that is two) number of plus indices. Distinguishing unpolarized (O), longitudinally polarized 
(L) and transversely polarized (T) situations (spin 1/2), 



-p+fj,;+:y -p+fj.: + u . -p+fj,:+u 



i rp , 



(Al) 



one has, with the (dimcnsionless) invariants a — 2k ■ P/M'^ and r = fc^/M^, the result 



-p+pi+i' 

^ o — 



P+ 

M 

P+_ 4^nr> 
M 



(-.9^) [(^2 + xX^) + x{Xi + XX3)] + 



P^ 



f^T '^n^ 



M 



{Xi + xX^ 



x]Xe, 



M J IvP 

P+ kl^n"^ 

iXf, 



Xe 



M M 



+n^'v}' 



r+^'+'' = iS,. 



2Xi + [(X2+xX4) + a;(X4 + .TX3)] + 2(|-a;)(X4 + xX3) + (|-x) Xg , 



2 ^krllJ. j^^f} 






-{Xis+X22)-{^-x]X2i 



M 



M 



X 
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--.X) {X15+XX13) 



(I - x) (Xi9 + X23) - (X21 + a;X23) - (I - a;) ' X25 



+5', 



M M 



n+/j;+i^ 



M 



M 



(I - x) (X18 + X22) - (^20 + a:X22) - (I - x)' 
r [-(^12 + a;Xi4) - x{Xi5 + xXia) 



X2 



'(Xi7 + xXiq) + (X21 + xXza) + ^ X25 






Af 



Af2 



+i 



J I t T' 't nj'T' ' i^ np 



M 



M 



M 



M M 

-(Xi5+a;Xi3) + Xi9 



^ ^ (X2o + a;X2 



StIa'i.''} 



Af 



(X16 + xXis) 



(f 



-a; X 



L25 



(A2) 



(A3) 



St [a* 1^] 
-^T n_ 



(Xg + xXii) + (I - x) (Xi7 + xXig) 



A//2 



X 



19 



1\7 






fC-p ' tJ 'jp ( ^ 



M 



(I - x) X24 + ef^^'^n!^ (I - x) (X16 + xXis) 



^n^kTSTlt^ jj'^} 



X18 
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M 



X 



18, 



(A4) 



We need the soft parts integrated over the momentum component fc~. Upon integration over dk~ we find 
that the functions appearing in the decomposition of the quantity MF*-' in Eq. 12 can be expressed in the 
amphtudes of Eq. |^. For this we use 



Mr^{x,kr;P,S) = Jdk- r+'^'+^k;P,S) = ^ f [dadr] r+''+^ {k;P, S), 



where we have introduced the shorthand notation 



[dcrdr] = dadr 6 [ t — xa + x 



AP 



(A5) 



(A6) 



to indicate integration over a and r. The resuhs for the twist two distributions in terms of the amphtudes Xi 
are the foUowing: 



[(^2 + xXi) + x{X^ + xXs)] + ^ Xe 



xG{x,k^) = / [dcrdr] 
vGrix, kl) = [[dadr] [{X^e + xX^s) - (^20 + .XX22)] , 



xH^{x,kl) = [dcrdr] Xq, 

xAGl{x, kl) = - [[dadr] [(Xg + xXia) + {Xq + xXn) + (| - x) [{X^ + xX^) + x{Xi5 + xXis)] 
kl 



xAGrix.ki 



M2 

[dadr] 



Xl9 +X23+{--x] X2S 



(X12 + xXii) + xiXi5 + xXu) - (Xn + xXig) - (X21 + xX^z) - j^ X25 



X AHtix, ki) = -2 / [dadr] (Xis + ^22) + ( - - a; ) X24 



xAH^{x,kl) = 2 [[dadr] X24, 
xAHT{x,kl) = 2 ([dadr] 



(X16 + xXis) + (X20 + XX22) - ^ X24 



From Eq. A2 we can also obtain the relation between the functions in the first quantity, MF' , relevant at 



twist three. We find that all functions have a real part (which is T-even) and an imaginary part (which is 
T-odd): 

Tie {xG^{x, kl)) = [[dadr] [(X4 + xX^) + (| - x) Xg] , 

Im{xGi:{x,kl)) = [[dcrdr] [-X5], 

Tie (xAGi^ix, kl)) = ([dadr] [Xio + (| ^ a^) (^i5 + xXi^) -{^-x) {X^ + ^23) 



-{X21 + xX2:i) - [- - x) X25 



- - a; ) (X18 +^22) + (X20 + .1X22; 



Im{xAG^^[x,kl)) = ([dadr] 

Ue (xAG3t(x, kl)) = I [dadr] [(Xg + xXn) + (f - 2:) {X^j + xX^^) 

'{Xi5 + xXii) - Xw 



2"^) ^24 



2M2 



2-2;)X25 
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Xm (^xAGsTixjk"^)) = [dadr] 



-(f-)(^ia + .X,,)-A_ 



X 



18 



7^e {xAG^Tix, kl)) = / [dadr] [-(X15 + xXia) + X19 + (^- - a; ) X25 
Xm(xAG^j,(x,fc^)) = /[dcrdT] 



--x]X24 



We now turn our attention to the different set of indices, namely (ij:l+), and the corresponding correlator 
r*-'^'+(fc; P, S). We again decompose according to the spin of the hadron and find 



-^ijd+ 



n^i,j] 



P+_ J -g'^'k 
Id 

P+ 



M 



[(X4 + xXs) - 1X5] 



r?''+ = ^ iS,el^ ^ 



M 



M 



Xll+{--x] (Xi4 + XX13) - (Xi7 + xX^g, 



x](X,g+X2 



(7 \ 2 



+S,e^ -^ 



M 



-(X16 + xXis) -{^-^) (^18 + ^22) - (I - 2; ) ^24 



(A7) 



-^ij;l + 



P+ 



le 



■ ■ n,.j. K'j' • (5^ 



M [ ' M M 



-{Xu + xXu) +X23+[--x) X25 



a 



f\j-p ni'j- ' iJ 'J- 



^ M M 



-X2 



x\X- 



(Xg + xXw) + (^1 - a;) (X21 + XX23) 



A^^'' 



el^SL 



(|-a;) (X2o+a;X22) + 



M2 



24 



X2 



Upon integration over dk , we find that the functions appearing in 



Mr^''(x,fcr)= / dk- r^''+(fc;P,5') 



2P+ 



[[dadr] r^''+{k;P,S) 



(A8) 



can be expressed in the amplitudes as follows: 

TZe{xHi-{x,kl)) = [[dadr] [X4 + XX3] 
Xm{xH^{x,kl)) ^ ([dadr] [-^5] , 



Tie {xAHi^ix, kl)) = I [dadr] [Xn + ( 2 - ^^j (^i4 + ^^13) - (^i7 + xX^g) 



-[--X^{X,g+X23)-[^-x) X25 



Im{xAH^^{x,kl)) = I [dadr] 



(X16 + xXi^) +['^-x) (X18 + X22) +('^-x] X2 



Ue [xAHtrix, kl)) = / [dadr] [-{X^a + xX^^i) + X23 + ( 2 " ""j ^25 

Im {xAHirix, kl)) = [[dadr] X22 - (| - 2;) X24 , 

7^e {xAHsxix, kl)) = Adadr] {(Xg + xXw) + (| - a;) (^21 + a;^23) 



+ 



2M2 



-(Xi4 + a;Xi3) -X23+{--x] X25 



Krj 



Xm {xAHMx, kl)) = J [dadr] \-[^- x) (^20 + XX22) - ^ 
The integrated functions in terms of the amplitudes are given by 



X- 



22 



2 - 2; ) ^24 



(A9) 
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xG{x) = M^TT / {dadr} 



{X2 + xXi) + x{Xi + xX^) + [xa -x^-t)^ 

xAGix) = -M^TT I {dadr} {(Xg + xXw) + (^9 + xX^) + {^ - x) [(X12 + xXu) + x(Xi5 + xXis)] 
-{xcr - x2 - r) Xig + X23 + (| - 2;) X25 } , 
7^e {xAG3t{x)) = M^n [{dadr} /.(Xg + xXn) + (^ - x) {Xn + ajXig) 

+ -(xcr - x2 - t) -(Xi5 + xXis) - Xig + (| - 2^) ^25 



Tee (a;Ai/3T(a;)) = M\ [{dadr} {(Xg + xXm) + {^ - ^) (^21 + a^Xsa) 

+ -(xcr-x2-T) -{Xii + xXu) - X23 + [^ - x^ X25 \, 
Im{xAH3T{x)) ^M^TT {dadr} I - (^^ - x'^ {X20 + XX22) - -(xcr - x^ - t) 



X 
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-x]X: 



i-24 



Xo 



x]X, 



with the convention 



{dcrdr} = dadr 9{xa — x — t). 



(AlO) 
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